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. Abstract 

. We present a new class of macroscopic models for pedestrian flows. Each individual 

is assumed to move towards a fixed target, deviating from the best path according to 
the instantaneous crowd distribution. The resulting equation is a conservation law with a 
nonlocal flux. Each equation in this class generates a Lipschitz semigroup of solutions and 
is stable with respect to the functions and parameters defining it. Moreover, key qualita- 
CN | tive properties such as the boundedness of the crowd density are proved. Specific models 

are presented and their qualitative properties are shown through numerical integrations. 
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1 Introduction 

From a macroscopic point of view, a moving crowd is described by its density p = p(t,x), so 
that for any subset A of the plane, the quantity J A p(t, x) dx is the total number of individ- 
uals in A at time t. In standard situations, the number of individuals is constant, so that 
conservation laws of the type dtp + div x (pv) =0 are the natural tool for the description of 
crowd dynamics. A key issue is the choice of the speed v, which should describe not only 
the target of the pedestrians and the modulus of their speed, but also their attitude to adapt 
their path choice to the crowd density they estimate to find along this path. 
Our starting point is the following Cauchy problem for the conservation law 



X: J d tP + dw[pv(p) {u(x)+l(p))j =0, (L1) 

c3 ■ \ P(0,x) = po{x) . 

The scalar function p i-> v(p) describes the modulus of the pedestrians' speed, independently 
from geometrical considerations. In other words, an individual at time t and position x G W N 
moves at the speed v (p(t, x)) that depends on the density p(t, x) evaluated at the same time 
t and position x. Given that the density is p, the vector v(x) +Z(p) describes the direction 
that the individual located at x follows and has norm (approximately) 1. More precisely, the 

individual at position x and time t is assumed to move in the direction v(x) + (z (x). 
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In situations like the evacuation of a closed space O C R , it is natural to assume that the 
first choice of each pedestrian is to follow a path optimal with respect to the visible geometry, 
for instance the geodesic. As soon as walls or obstacles are relevant, it is necessary to take into 
consideration the discomfort felt by pedestrians walking along walls or too near to corners, 
see for instance |19|. 121] and the references therein. 

The vector (x) describes the deviation from the direction v(x) due to the density 

distribution p(t) at time t. Hence, the operator X is in general nonlocal, so that {x (/?(£)) ^ (x) 
depends on all the values of the density p(t) at time t in a neighborhood of x. More formally, 
it depends on all the function p(t) E L 1 (M Ar ; [0, R]) and not only on the value p{t,x) € [0, R]. 
The case in which I = is equivalent to assume that the paths followed by the individuals 
are chosen a priori, independently from the dynamics of the crowd. 

Here we present two specific choices that fit in (jl.ip . A first criterion assumes that each 
individual aims at avoiding high crowd densities. Fix a mollifier r]. Then, the convolution 
(p * rf) is an average of the crowd density around x. This leads to the natural choice 

MP) = ~g l == , (1-2) 

'1+ ||V(/9*7/)f 

related to [3], which states that individuals deviate from the optimal path trying to avoid 
entering regions with higher densities. Through numerical integrations, below we provide 
examples of solutions to (|l.ip - (|1.2p . They show the interesting phenomenon of pattern for- 
mation. In the case of a crowd walking along a corridor, coherently with the experimental 
observation described in the literature, see for instance [T71 [TBI [201 [28] , the solution to 
(|1.2p self-organizes into lanes. The width of these lanes depends on the size of the support of 
the averaging kernel rj. This feature is stable with respect to strong variations in the initial 
datum and also in the geometry. Indeed, we have lane formation also in the case of the evac- 
uation of a room, when the crowd density sharply increases in front of the door. Section 14.11 
is devoted to this property. 

From the analytical point of view, we note that the convolution term in (|1.2p seems 
not sufficient to regularize solutions. Indeed, the present analytical framework is devised to 
consider solutions in L 1 RBV. Both in the case of a crush in front of an exit (Section I4.2p and 
in the specific example in Section |4. 31 numerical simulations highlight that the space gradient 
of p increases dramatically. 

According to (jl.2p . pedestrians evaluate the crowd density all around their position. When 
restrictions on the angle of vision are relevant, the following choice is reasonable: 

x(p)= c v I^p(y) r i( x -y) ( p{(y- x )-9ix))dy ^ 



i + 



V f RN p{y) r](x - y) tp((y - x) ■ g(x)) dy 



Here, 77 is a fixed mollifier as above and the smooth function <p weights the deviation from 
the preferred direction g(x). 

We note that the constructions in [221 123j and [T4J fit in the present setting. Indeed, there 
the following choices were considered: 

/ dt P - div (p f (p) Vtp) = / d t p - div (p f 2 (p) Vip) = 

1 |V^| = l//(p) and 1 - e A^+|V^| 2 = l/(/(p) + e) 2 . 
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The former admits an immediate interpretation: the direction V(p of the speed v is chosen 
along the solutions to the eikonal equations, i.e. all pedestrians follow the shortest path, 
weighing at every instant the length of paths with the amount of people that are moving 
along it. In the former case, v{x) = and 1{p) is the gradient of the solution to the eikonal 
equation with l/f(p) in the right hand side, while in the latter case v(x) = and I(p) is the 
gradient of the solution to an elliptic partial differential equation. 

The model introduced in [13} [29] relies on this measure valued conservation law: 

dtp + div (p v) = where v = v{x) + I f(\x — y\) (p ((y — x) ■ fix)) j rrdp+iy). 

Jr 2 If- 2/11 

Here the unknown is a map p : [0, T] — > A4(W N ; M) assigning at every time t a positive measure 
p(t) which substitutes the crowd density, in the sense that the amount of people that at time 
t are in A is (/^(t)) (^4). Contrary to this model, the present framework is the space of L 1 
densities and physical a priori L°° bounds on the solutions to (|l.ip are rigorously proved, 
preventing any focusing effect as well as the rise of any Dirac delta. Moreover, below we prove 
global in time existence of solutions, their continuous dependence from the initial data and 
their stability with respect to variations in the speed law v. 

In |12| . the geometric part v of the speed is chosen a priori, while its modulus depend on 
the density as well as on the gradient of the density: 

d t p + div (p (p(p, Vp) v{x)) = . 

On the contrary, the model introduced in [7] postulates a nonlocal dependence of the speed 
from the density: 

d t p + div (pv(p * rf) v{x)) = , 

which amounts to assume that pedestrians choose their behavior according to evaluations of 
an average of the density around their position, rather than according to the density at their 
place. The following second order model was presented in [HUB] and does not fit in (j 1 . 1 p : 

dtp + div (pv) = 
d t v + (v ■ V)v = F(p, Vp, v) , 

we refer to the review |2J for further details. To underline the variety of analytical techniques 
with which crowd dynamics has been tackled, we recall the further approaches: optimal 
transport in [S], the mean field limit in [IS], the functional analytic one in [26, .27] and the 
nonclassical shocks used in [6j [TUl E] ■ 

Throughout this work, we set all statements in all R . In Section [3] we show how the 
present framework is able to take into consideration the presence of various constraints. In 
the case of pedestrian dynamics, for instance, this amounts to prove that no individual passes 
through the walls of a given room, provided the initial datum is assigned inside it. Concerning 
the dimension, our main applications are referred to the case of crowd dynamics, i.e. N = 2. 
Nevertheless, from the analytical point of view, considering the case of a general N does not 
add any difficulty. Furthermore, we believe that the present setting can be reasonably applied 
also to, say, fishes and birds moving in M 3 , for example in predator-prey like situations as 
described in |8J. 

The next section is devoted to the analytical properties of well posedness and 

stability. The general theory is particularized to specific examples in Section [3] where the 
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presence of obstacle (walls) is considered. Sample numerical integrations are provided in 
Section UJ The final Section collects the analytical proofs. We defer to Appendix lAl further 
remarks related to the geometry of the physical domain. 



2 Analytical Results 

In the following, N G N \ {0} is the (fixed) space dimension. We denote R + = [0, +oo[; 
the open ball in centered at x and with radius r > is B(x,r) and we let Wn = 

f* /2 (cos #) N dtf. As usual, we denote S^" 1 = {x G R N : \\x\\ = l}. 

The density p can be defined as the fraction of space occupied by pedestrians, so that p 
turns out to be a nondimensional scalar in [0, 1]. Otherwise, it can be useful to think at p 
as measured in, say, individuals/??! 2 and varying in [0, R], with R > being a given maximal 
density, for example R = 8 individuals/m 2 . 

Our first step in the study of (jl.ip is the formal definition of solution. 

Definition 2.1. Fix a positive T and an initial datum p$ G L 1 (IR Ar ; [0, R]). A function 
p G C° ^[0, T]; L 1 (M Ar ; M)^ is a weak entropy solution to (jl.ip if it is a Kruzkov solution to 
the Cauchy problem for the scalar conservation law 

J dtp + div (pv(p)w(t,x)) =0 . . 

\ p(0,x)=po(x) 

where w(t,x) = u{x) + (z(p(i))^ (x). 

In other words, recalling [Ml Definition 1], for all k € R, for all ip £ C^°(]-oo,T] x R iV ;M + ), 

n\p — k\ dt(f + (pv(p) — k v(k)) w(t, x) ■ V<p — k v(k) div w(t, x) ip sgn(p — k) dxdt 

+ / \pq(x) — k\ ip(0,x) dx ^ 0. 
Jr n 

On the functions defining the general model (jl.ip . we introduce the following hypotheses: 
(v) v G C 2 (R; R) is non increasing, v(0) = V and v(R) = for fixed V, R > 0. 
(i/) G (C 2 n W 1 '°°)(R Ar ;R iV ) is such that div v G (W 1 ' 1 n W 1 '°°)(R iV ; R). 

(I) X G C° (VqR^; [0,i?]);C 2 (R iV ;R iV )) satisfies the following estimates: 

(1.1) There exists an increasing C/ G L^ C (R+;R+) such that, for all r G L^R^; [0,R]), 

||^)|| wl ,oc < Cj(||r|| Ll ) and ||divJ(r)|| Ll < Cj(||r|| Ll ) . 

(1.2) There exists an increasing d G L^ C (R + ;M + ) such that, for all r G L 1 (R JV ; [0,R]), 

||VdivX(r)|| Ll <C / (||r|| L1 ). 
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(1.3) There exists a constant Kj such that for all r\,r 2 G L 1 (M Ar ; [0,/?]), 

||2:(ri)-I(r 2 )|| L0O ^ K/-|ki-r 2 |I Ll , 
\\X{r l )-X{r 2 )\\^ + \\diY{X{r 1 )-X{r 2 ))\\^ s$ Kj ■ ||n - r 2 || L i . 

Furthermore, throughout we denote by q the map q: [0,R] t— > R defined by g(p) = pv(p). 
As a first justification of these conditions, note that they make Definition 12, II acceptable. 

Lemma 2.2. Fix a positive T. Let (v), (v) and (1.1) hold. Choose an arbitrary density 
r G C° ([O.TJiL 1 ^; [0,i2])V Then, the Cauchy problem 



{ X) d = & w(t,x)) = mtk w{t ' x) = v{x) + ( i(rW) ) {x) (2 - 2) 

satisfies the assumptions of Kruzkov Theorem \24[ Theorem 5]. 

This lemma is proved in Section [5] below. In Section [3] we show that the above assump- 
tion (I) allows to comprehend physically reasonable cases. 

The next result is devoted to the proof of existence and uniqueness for (jl.ip . It is obtained 
through Banach Fixed Point Theorem. 

Theorem 2.3. Let (v), (u) and (I) hold. Choose any p G (L 1 nBV)(R iV ; [0,R]). Then, 
there exists a unique weak entropy solution p G C° ^R + ;L 1 (R Ar ; [0, R])\ to Moreover, 
p satisfies the bounds 



||p(*)|| L i = IIPo|| L i . f° r a - e - 1 e 

'([0.JJJ) 

where k = (2N + l)||g'|| LOO([0ijR]) (||Vi/|| Loo + Cj(||po|| l i)) 



TV(p(t)) < W(p )e w + te w W Ar ||g|| LooM) (||Vdiv l /|| Ll +C'i(||po|| L i)) , 



Using the techniques in [71 [9] , we now obtain the continuous dependence of the solution 
to (jl.ip from the initial datum and its stability with respect to v, v and X in the natural 
norms. 

Theorem 2.4. Let (v), (v) and (I) be satisfied by both systems 

dtp + div pvi{p) (vi(x) +X l (p)) =0 Idtp + div pv 2 (p) {v 2 {x) +X 2 (p)) 

p(0, x) = p ,i(x) |p(0, x) = p ,2(x) 

with po,i,pa,2 G (L 1 nBV)(R w ; [0, R]). Then, the two solutions p\ and p 2 satisfy 

\\Pi{t) - P2{t)\\- Ll ^ C(t)(\\p 0jl -p ,2|| L i + lift -ffellw 1 . 00 

+ 11^1 - ^2|| L oo + ||div(i>i - ^2)|| L i +d(Xi,X 2 ) 

where d(X x ,X 2 ) = sup {\\li(p) - Mp)\\^«> + ||divXi(p) - divZ 2 (p)|| Ll : p G L 1 ^; [0, 12])}, 

the map C G C°(R + ;R + ) vanishes at t = and depends on TV(poi)> po,i Ld II^iIIl 00 ' 
||div 

^lllw 1 ' 1 ' lift II W 1 ' 00 ' 11^2 Hw 1 ' 00 ' 
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Thanks to these stability results, several control problems can be proved to admit a 
solution through a direct application of Weierstrafi Theorem. A possible standard application 
could be the minimization of the evacuation time from a given room. Describing the actions 
of a controller able to determine the initial pedestrians' distribution and/or their preferred 
paths, one is lead to an optimal control problem with the initial datum po and the vector 
field v as control parameters, for instance. Without any loss of generality, the compact sets 
on which the optimization is made can be |/9 € (L 1 n BV)(lR iV ; [0, 1]) : TV (p) < Af| for po 

and ji/ G C 2 (r n ; B(0, 1)) : ||V 3 z^|| < m\ for u, where M > is arbitrary. 

A different problem solved by the same analytical techniques is that of the dynamic 
parameter estimation. Once real data are available, one is left with the problem of determining 
the various parameters entering v, v or X. Theorem [23] ensures the existence of the parameters 
that allows a best agreement between the solutions to (|l.ip and the data. 

3 The Models 

This section is devoted to the study of specific cases of (|l.ip . Aiming at real applications, it is 
necessary to take into consideration the various constraints that are present on the movement 
of pedestrians. Therefore, we introduce the subset Q, of which characterizes the region 
reachable to any individual. The boundary d£l consists of walls that can not be crossed by 
any individual. On the set f2 we require that 

(ft) C R w is the closure of a non empty connected open set. If d£l is not empty, there 
exists a positive such that the function 

d dn : B(dn,r Q )nn -> r+ 

x i-> inf {d(x, w) : w € <9$7} 

is of class C 2 (B(dn, r n ) n O; M+) . 

Note that we do not require £1 to be bounded. This assumptions allows us to introduce the 
inward normal 

n(x) = Vd dQ (x) (3.1) 

on all the strip £?(<9f2,rn) n Q. Moreover, ||n(x)|| = 1. 

For the present models to be acceptable, it is mandatory that no individual enters any wall, 
provided the initial datum is supported inside the physically admissible space. Analytically, 
this is described by the following invariance property: 

(P) The model (jl.ip is invariant with respect to SI if 

spt po C fl spt p(t) C f2 for all t ^ . 

Below, we show that theorems 12.31 and 12.41 can be applied to reasonable crowd dynamics 
models, in the sense that the presence of a physically admissible set f2 is also considered and 
its invariance in the sense of (P) is proved. More precisely, we show below the following 
sufficient condition for invariance. 
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Proposition 3.1. Let u, X and satisfy (u), (I) and (il). If for all x € d£l and p € 

L 1 (R 7V ; [0,R]) with sptp C 

(y(x) + (X(p)) (s)) • n(x) > (3.2) 

then, property (P) holds. 

3.1 The Model f lLTf - fO]) 

The starting point for (|l.ip - (jl.2p is provided by the degenerate parabolic model 

/ / 



dtp + div 



pv(p) \ u(x) - 



y \ ^/i + ||vv(p) 
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introduced in |4j, which fits in (|l.ip with I(p) = —eVip(p) I \J f + || V^(p)|| 2 , motivated 

by the desire of each individual to avoid entering regions occupied by a high crowd density. 
Here, ip is a suitable weight function. Assuming that each pedestrian reacts to evaluations of 
averages of the density, we obtain 



d t p + div 



/ 

mil n\ \ u( nr\ — &■ 

\ - ^/l + HvtiKp)*!)) 




Here, we avoid the introduction of ip to limit the analytical technicalities, obtaining (jl.ip - 
(|1.2p . We assume throughout that e > is fixed and that the mollifier n satisfies 

(rj) Tj G C3(R JV ;]R+) with j RN r)(x) dx = 1. 

This convolution kernel has a key role: the individual at x deviates from the optimal path 
considering the crowd present in the region x — spin, when no walls are present. The value 
r](£) is the relevance that the individual at x gives to the density p(x — £) located at x — £. 

The next result shows that the present model fits in the framework described in Section [2] 
when applied on all of W N . 

Lemma 3.2. Fix e > and let rj satisfy (r\). Then, the operator X in (jl.2p satisfies (I). 

When the region 17 reachable by the crowd is restricted by the walls d£l, we intend the 
convolution restricted to f2 

{p*v)(x)= / p(y) v(v - x ) d v ( 3 - 3 ) 

Jn 

which coincide with the previous definition in the case spt p C 0. A better choice is described 
in (jA.ip . The vector v is here chosen as a sum 

u = g + S. (3.4) 

The former vector g is tangent to the "optimal" path, depending only on the geometry of the 
environment and coherent with it. Hence we assume that 
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(g) g G C 2 (R N ;S N ~ 1 ) satisfies Vg E L 00 ^; M Arx7V ), dive/ G (W 1 ' 1 n W 1 - 00 )^; M) and 
the invariance condition g(x) ■ n(x) ^ holds for all x G d£l. 

The latter vector 5 describes the discomfort felt when passing too near to walls or obstacles. 
Below, we choose 

S(x) = A a(x) n(x) . (3-5) 

Here, A G M + is a suitable constant and n(x) is the inward normal (|3.ip . The role of the 
function a is to confine this discomfort to the region near the walls, i.e. we require 

(a) a G C 2 (M Ar ; [0, 1]) is such that a(x) = whenever B(x, rpj C and a(x) = 1 for 

x G da 

A better choice for the discomfort is discussed in Appendix [Al 

The present setting (|l.ip - ([1.2p can be effectively applied also in presence of walls, as 
shown by the next result. 

Proposition 3.3. Let e > 0, (v), (SI), (rj), (g) and (a.) hold. Define v by (|3.4p . 5 
by (|3.5|) and Z by (|1.2p . Assume moreover that either d£l is compact, or a G (W 2 ' 1 n W 2,oc ) 
and don G (W 3 ' 1 n W 3 '°°). Then, (|l.ip - (jl.2p satisfies the assumptions of Theorem [2731 and 
Theorem \2.4\ Furthermore, if A ^ ei? ||Vr/|j L i ; then property (P) holds. 

3.2 The Model f lLT|) - ffl~3]) 

In the framework of (jl.ip . we now extend the model |13|. (2.1)— (2.4)] to take into account 
the effects of crowd density on the speed modulus, i.e. of v(p). Moreover, we interpret the 
nonlocal term as a weighted average of the density, where the amount of crowd in the direction 
of the optimal path g is given more importance. 

Using the same notation as above, we thus consider the model (jl.ip with 

_ _ £ VJrn p(y) y(x - y) <p((y - x) ■ g{x)) dy ^ 

1+ v lRNp(y)v( x -y) p[(y-x) ■ g{x)) dy 

Here, e, n and v are as in § 13.11 and, in particular, v is as in (|3.4p . 

The (almost) isotropic convolution (|3.3|) is here weighted by <p, whose argument is essen- 
tially the angle between the preferred path g and y — x, the point x being the position of 
the individual an the point y being the location where the density is evaluated. This term 
takes into account the preference of each individual to deviate little from the preferred path 
defined by g. For example, let (p G C°°(1R, [0, 1]) be such that ip = on ]— oo, 0] and tp = 1 on 
[■#, +oo[, where $ > is a given parameter. Then, adding the function tp ((y — x) ■ g(x)) into 
the nonlocal term means that the individual at x reacts to the average density evaluated in 
the prefer ed direction g(x). 

The denominator is a regularized normalization. Its presence is necessary from the model- 
ing point of view, for coherence with the presence of v(p). From the analytical point of view, 
this normalization makes various expressions slightly more complicate, but all estimates re- 
main doable. 



S 



We are thus lead to the equation 
/ / 

V/ Rjv p(y) r)(x -y)ip((y-x)- g(x)) Ay 



d t p + div 



pv(p) 



\ 



v{x) 



\ 



1+ ^Irn p(y)v(x -y)f((y - x) ■ g(x)) dy 



(3.7) 



n 



and we verify that it fits in the analytical framework provided in Section [2j 

Lemma 3.4. Fix e > 0, 77 G C;!(R JV ;M+) with J RN r](x)dz = 1, g G W 3 >°°(R; [0, 1]) and 
v G W^R^R). Then, the operator 1 defined by JOJ satisfies (I). 

Passing now to the case in which the crowd's movement is constrained by the walls <9f2,as 
in the preceding section, we intend all integrals in (|3,6p - (|3.7p restricted to O. In particular, 
we consider now 



ZOO 



Y/o p(v) v( x - y) ^{(y - ■ g{ x )) dy 



(3.8) 



1 + 



V Jo P(v) V(x - y) <p[{y - x) ■ g{x)) dy 



with v as in (v) and v as in (|3.4p . The applicability of theorems 12.31 and 12.41 to (|l.ip - (|3.8p 
and the validity of property (P) is ensured by the following proposition. 

Proposition 3.5. Let e > 0, (v), (SI), (17) and (g) hold, with moreover g G (C 3 f~l 
W 3 ' 0O )(M JV ;5 JV - 1 ). Let^e (C 3 n W 3 >°°)(R; R). Define v by d33} 3 6 by ([32]) and 1 by f^Hjl . 
Then, (jl.ip f)3.6|) satisfies the assumptions of Theorem \2.3\ and Theorem \2.4\ 
Moreover, assume that ip' and call i = diam(sptry). Then, 



ip' and 
imply that property (P) holds. 



A ^ -R IMIw 1 ' 1 IMIw 1 ^ (1 + ^ ||Vfl|| Loo ) 



3.3 The Model [13], (2.1)-(2.4)] 

The model in |13|. I29j. although set therein in the space M(WL N ;~. 
sures, can be seen as a particular case of (jl.ip setting 



") of positive Radon mea- 



v(p) 
v(x) 



1 

^des (x) 



and l(p) =£ p(y) Vrj(x - y) ip((y - x) ■ u(x)) dy , (3.9) 



v das being the desired speed, see [131 formula (2.4)]. When rj is radially symmetric, we re- 
cover [131 formula (2.6)] with rj(x) = 7?(||x||) and fj' = f. We leave to [13] the motivations of 
this model. 



Proposition 3.6. Let e > 0, (v), (v) and (rj) hold. Assume that <p G (C 2 n W 2 '°°) 
Then, (jl.ip ~ (|3.9p satisfies the assumptions of Theorem \2.0A and Theorem \2.4\ 

In particular, above we prove that if the initial datum po is m L 1 (R Ar ; [0, it!]), then the 
corresponding solution satisfies the same bounds. This ensures that neither focusing to any 
Dirac delta takes place, nor values of the density above R can be expected. (For the sake of 
completeness, we note that the conditions r\ ^ and J RJV rjdx = 1 is in the case of (jl.ip — f|3.9[) 
neither necessary, nor meaningful and can be replaced by r\ G (C 3 n W 2, °° n W 3,1 )(R 7V ; R), 
see the proof in Section [5] for more details.). 
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4 Qualitative Properties 



This section is devoted to sample numerical integrations of In all the examples be- 

low, we choose as vector field v = v(x) the geodesic one, computed solving numerically the 
eikonal equation. This leads to the formation of congested queues near to the door jambs. 
Prom the modeling point of view, a more refined choice would consist in choosing v so that 
most pedestrians are directed towards the central part of the exit. This choice increases the 
difficulties neither of the analytical treatment nor of the numerical integration but imposes 
the introduction of several further parameters. 

The algorithm used is the classical Lax-Priedrichs method with dimensional splitting. As 
usual, a fixed grid (xi,yj) for iq, . . . ,n x and j = l,...,n y is introduced and the density p 
is approximated through the values pij on this grid. At every time step, the convolution in 
vector I(p) is then computed through products of the type A^p^kB^j, for two fixed matrices 
A and B depending only on the kernel r\. 

All the examples below are set in M 2 , due to obvious visualization problems in higher 
dimensions. As is well known, the analytical techniques are essentially independent from the 
dimension as also the numerical algorithm. The time of integration obviously increases with 
the dimension. 

4.1 Lane Formation 

A widely detected pattern formed in the context of crowd dynamics is that of lane formation, 
see for instance [T71 [THJ [201 12H] • This feature has been often related to the specific qualities 
of each individual, i.e. it has usually been explained from a microscopic point of view. Here, 
in a purely macroscopic setting, we show that the solutions to (|l.ip - (|1.2p also display this 
pattern formation phenomenon, with pedestrian self organizing along lanes. 



Time- 0.000 



Time- 2.529 



Time = 5.043 







Time = 7.557 



y 



Time- 10.071 




D 



Time - 15.014 



D 



Figure 1: Solution to l^-(^~<^~§FQ) at times t = 0, 2.529, 5.043, 7.557, 10.071, 15.014. 
First 3 lanes are formed, then the middle lane bifurcates forming the fourth lane. 
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Consider (frTj) - (fT2]) with 



v{x) 
v{p) 



+ 8{x) , rj(x) 



\ (1 - p) , p (x) 



X 



l-m 2 ] [l-(f) 2 ] 
x) , 



[3/5,4] x [-3/5,3/5] 



4 

5 ' 



(4.1) 



where 5 = 5(x) describes the discomfort due to walls: it is a vector normal to the walls, 
pointing inward, with intensity 3/2 along the walls, decreasing linearly to at a distance 
3/10 from the walls. As Figure [T] shows, the initially uniform crowd distribution evolves into 
a patterned configuration, first with 4 lanes and then with 5 lanes. The number of lanes 
depends on the size of the support of the convolution kernel rj. Indeed, keeping all functions 
and parameters fixed, but not the parameter r, which determines the size of sptr/, we obtain 
patterns differing in the number of lanes, see Figure 



1 = 9.014 1 = 9.014 1 = 9.014 






D 




Figure 2: Solution to (fT^ ~ (fL2]) - (^I]) - (^T]) computed at time t = 9.014 and with sptr/ with 
radius r = 0.5, 0.6, 0.8, 0.9, 1.0, 1.4. Note that as r increases, the number of lanes diminishes. 

The formation of lanes is a rather stable phenomenon. Indeed, Figure [3] shows the result of 
the integration of <^2)-<^2i)- §53$) -dfCU computed at time t = 0, 5, 10 with r = 3/5 (above) 
and r = 9/10 (below) with initial data different from that in (j4. 1[) . In both cases, lanes are 
formed similar to the corresponding situations in Figure [2j 

We also note that in the present framework, using the terms in [281 Section 5.4], lanes 
form also in an isotropic setting. Indeed, the integrations in figures HH2] were obtained with 
individuals able to see both forward and behind. 



4.2 Evacuation of a Room 

A standard application of macroscopic models for crowd dynamics is the minimization of 
evacuation times. The present setting applies to general geometries, see the assumption (f2). 
Here we show that (jl.ip captures reasonable features of the escape dynamics. 
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D 




D 



D 



3 





D 



Figure 3: Solution to (fT3]) -(fT ^ - (|3T4^ - (j4^ with different initial data at time t = 
0, 5.014, 10.014 and above with r = 0.6, below with r = 0.9. Note that above 4 lanes 
form and below 5, similarly to what obtained in Figure [TJ 



We consider a room with an exit, as in Figure [U The vector v = v{x) is chosen as the 
unit vector tangent at x to the geodesic connecting x to the exit. The discomfort d = d{x) 
is a vector normal to the walls, pointing inward, with intensity 1 along the walls, decreasing 
linearly to at a distance 1/2 from the walls. The other quantities are in (|4.2p . Keeping the 



Time= 0.000 Inside = 100.00% 



| 
I 



v{p) 


= 6(1 -p) 






Tj(x) 






3 




r 


= 0.6 






Po{x) 


= 0.75 x, , 

[2,7] x 


-2,2] v ' 


e 


= 0.4 







(4.2) 



Figure 4: Initial datum and room 
geometry considered in § 14.21 

above parameters fixed, as well as the outer walls of the room, we insert various obstacles 
(columns) to direct the movement of the crowd. 

First, Figure \5\ first line, shows an integration of the case with two columns that direct 
people towards the exit. The number of lanes self adapts to the available space, with three 
lanes merging into one before the bottleneck. On the second line of Figure [5J the insertion of 
three columns in these positions delays but does not avoid the congestion at the exit. Those 
individuals that pass through the bottleneck are favored in exiting the room. 

In Figure El first line, the insertion of four columns is more successful. Note, in the first 
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Time = 2.521 Inside = 75.20% 



Time = 5.043 Inside = 46.85% 



Time- 7.563 Inside = 21.41% 




Figure 5: Solution to (|l.ip - (|1.2p - (|3.4p - (|4.2p with different geometries, computed at time 
t = 2.521, 5.043 and 7.563. 



two diagrams, that the number of lanes changes from 7 before the bottleneck, to 4 in between 
it. The individuals that do not pass through the bottleneck are penalized by the high density 
they meet near the door jambs. Finally, on the second line in Figure El an asymmetric layout 
hinders the lanes' pattern. 

4.3 On the Rise of Singularities 

A typical feature of conservation laws is the possible rise of singularities, see for instance [3j 
Example 1.4]. The nonlocal equation (jl.ip shares this characteristic. Indeed, assume p = 
p(t,x) is a given solution to (jl.ip . smooth up to time T > 0. Then, setting w(t,x) = 

(v(x) + (l(t)) (x)j , simple computations lead to the following equation for the space deriva- 
tive pj of p in the direction Xj, for j = 1, . . . , N: 

d t pj + (div pj) q w = (pj) 2 q" w + pj q" '^(dip) w + pj q div w + q '^{dip) djiv + q dj div w . 

The first term in the right hand side is quadratic in pj, showing that a blow up of pj may 
take place in finite time. 

Moreover, we consider (jl.ip with I as in (jl.2p and 

v(p) = 1-p, v(x,y) = (1 - 16a; 2 ) 3 (l - 16y 2 ) 3 x [1/41/4]2 (x,y), 

u(x) = 0, Po (x,y) = (l-4x 2 /9) 2 (l-4y 2 /9) 2 X[ _3 A 3 /2]2 (x,y), £=_ ' 

(4.3) 

and we obtain the solution in Figure [7] 
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Time = 2.521 Inside = 76.62% 



Time = 5.043 Inside = 48.64% 



Time = 7.563 Inside ^ 23.64% 





U 



Time . 2.521 Inside . 75.82% 



Time = 5.043 Inside = 47.36% 



^3 fl 



Time- 7.563 Inside = 21.85% 




Figure 6: Solution to (jTI|) -(fP j) -(|g^ j) - 
t = 2.521, 5.043 and 7.563. 



with different geometries, computed at time 




Figure 7: Left, the initial datum and, right, the solution to (|l.ip - (|1.2p - (|4.3p at time t = 1. 
Note the formation of vertical faces in the initially smooth distribution. 
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5 Technical Details 



5.1 Existence and uniqueness 

Proof of Lemma \2.2l Let q(p) = pv(p). Thanks to (u), the assumptions on r and (1.1), 

we have for all M > 

d p f{t,x,p) = q'(p) (v(x)+l(r(t))(x)); d p f G L°°([0, T] x R N x [-M, M]; R N ) ; 
dwf(t,x,p) = q(p) div (u(x)+l(r(t)) (a?)) ; div / G L°°([0,T] x I N x [-M,M];R). 



Thus, we can apply Kruzkov Theorem \2A\ Theorem 5 & 5.4] and ensure that (|2,2p admits 
a unique Kruzkov solution p G L°° (r + ; L 1 (R n ; R)J , which is continuous from the right in 
time. □ 



Lemma 5.1. Let p G L 1 ^; [0,12]), r G C° (m+jL 1 ^; [0,i2])j. Under hypotheses (v)- 
(u)-(I.l), the Cauchy problem (|2.2p admits a unique weak entropy solution p, with p G 
C° (R+; L 1 {R n ; [0, R])\ satisfying, for all t G R + , 

||p(*)|| l i = IIPo|| L i • 

Assume, in addition, that (1.2) is satisfied. Then, po G BV(R Ar ; [0, 12]) implies p{t) G 
BV(R' V ; [0, R]) for all time t ^ 0. Moreover, the following bound is satisfied 

TV(p(*))< TV(p )+tW^||g|| Loc(M) (||Vdiv l ,|| L1 +C / (||r|| Loo([0iT];L1) )) e^, (5.2) 
where Wn = J^ 2 (cos , d) N di? anc? i/ie constant Kq is bounded above as follows: 
k* < (2A r + l)||g'|| Loo([0jR]) (||Vi/|| LOO +C/(lkll L -([o,T];Li))) • 

Proof. Below we denote g(p) = pv(p), so that f(t,x,p) = q(p) (v(x) + (l(r(t)^ {x)\ . 

The existence of a solution follows from Lemma [2.21 The rest of the proof is obtained through 
the following steps. 

Estimates in L 1 and L°°. Since p is a solution to a conservation law, (|5.ip is immediately 
satisfied. 

Besides, p = and p = R are solutions to (|2.2p associated respectively to the constant 
initial conditions po = and po — R- Hence, we can apply the comparison Theorem [24, 
Theorem 3] and obtain that if the initial condition takes value in [0,12], then the solution 

takes values in [0, 12], so that p G L°° (r + - L 1 ^; [0, R]f 



15 



Continuity in time. Thanks to [9j Remark 2.4] or [25, Corollary 2.4], p is continuous in 
time if, for any T > 0, 



MIl°°([o,.r]) 



'o Jr n 

We have, thanks to (1.1), for all t ^ 
div (v+l(r(t) 



div z/(x) + (l(r(t)) (x 



dx dt < oo . 



(5.3) 



< l|divf|| L i+ C/(||r(t)|| Ll ). 

L 1 

Hence, condition (|5.3p is satisfied under hypotheses (v)-(i/)-(I.l) and we also obtain p G 
C°(R+;L 1 (K JV ;[0,i?])). 

Estimate in BV. To prove the bound on the TV norm we use [25} Theorem 2.2]. To this 
aim, we have to check that for any T > 0, 



Vd p f eL°°([0,T}xR N x[0,R};R NxN ) and f [ ||Vdiv/(i,x, 

Jo Jr n 



lL°°([0,.R]) 



dx dt < oo . 



Note first that Vd p f(t,x,p) = </(/?) V [v + X(r(t)) ). Hence, thanks to (v)-(i/) and (1.1), 
we have 



\Vd p f(t,x,p) 



IL°° 



Il°°([o,-R]) 



|Vi/|| LOQ + C7j(||r(t)|| Ll )) , 



Thanks to (v)-(i/) and (1.2), we have, for all i ^ 



Vdiv(i/ + Z(r(t))) < ||Vdivi/|| L i + Cj(||r(t)|| Ll ) . 
Applying |25(. Theorem 2.2], we obtain 

k* = (2iv + i)||va p /|| LOO 

< ( 2N + 1) Ik'IL-ao.il]) (H W Hl~ + C/(lkllL»([0,3l;Li))) > 

TV(p(t)) < TV(p )e K o t 

+te K ° t NW N \\q\\ Loc{m]) (||Vdivi/|| L i + C / (||r|| Lo o ([0)T] . L i ) ) 

completing the proof. 



□ 



Proof of Theorem \2.3l Thanks to Lemma 15.11 under hypotheses (v)-(iv) and (I), for 

any r G C° ([0, T];L 1 (R JV '; [0, i?])) there exists a unique solution p G C° (lR+; L 1 (R iV ; [0, i?])) 
to (|2.2p satisfying (|5.ip - (|5.2p . For any time T > 0, introduce the space 

X(T) = jr G C° ([O^L 1 ^; [0,fl])) : Vt G [0,T], ||r(t)|| L1 < ||p || L i} • 

Note that X(T) is a Banach space with respect to the norm |MIx(t) = su Pte[o,T] || r (*)||L 1- 
Let Q be the map that associates to r G X(T) the solution p G X(T) to (|2.2p . We want to 
find a condition on T so that Q is a contraction. 
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Let ri,r 2 £ X(T). To apply |25^ Theorem 2.6], we have to check that, for all t E [0, T], 



d P {h ~ h) E L°°([0,T] xi^x [0,i?];lO , and 



|div(/i-/ 2 )(t,a;,-)| 



L°°([0,fl]) 



dxdt < +00 , 



Thanks to (1.3) we have 



9 P (/i-/a)|| Loo = (X(ri)-Z(r 2 ))g' < K T ||g'|| Loo{[0)fl]) ||n - r 2 || Loo([0iT];Ll) 



|div(/i-/ 2 )(*,x,.)| 



/0 JM^ 

Hence, we get for all t ^ 



L°°([o,fl]) 



dxdt ^ TKj |[g|| L °°([o,R]) l! r i ~~ ^IIl^qo.tIjL 1 ) 



||(pi -p2)(i)|| L i 

^ * ^lkllw 1 >°°([0,i?])li r l ~ ^^^([O.TJjL 1 ) 



1 + e K ° ( TV (po) H —— t lkllL°°([o,_R]) 



|Vdivi/|| L i + C/(||r|| LOO([0)T] . L i)) 



Denoting 



|vhIl«+C/(II«)IIlO) 



(5.4) 



* = (2A^ + i)||q'|| loo([0R]) 
c 2 = -^-^||g|| L -([o,fl]) (ll vdivi/ llLi + C/(IIpo|| l i)) • 

The above estimate can be written 

||(Pi-p2)(i)|| Ll «S Cxt\\r x -r 2 || LOC([0)T];Ll) (l + e kt (TV (p ) + C 2 t 
We choose now T so that 

d r (l + e fcT (TV (p ) + C 2 T)) = 1 

and, applying Banach Fixed Point Theorem, we obtain a unique fixed point for Q on X(T). 
We now iterate the procedure above. To this aim, we use the total variation estimate (|5.2p 

TV(p(i)) ^ (TV(p ) + 2C 2 t)e kt , 

and, given T n , we recursively define T n+ \ so that 

Ci (T n+1 - T n ) (l + e^ 1 TV (p ) + C 2 e fe ( T «+i- T «)(2T n e fcT « + T n+1 - T n 



and the above procedure ensures the existence of a fixed point on the interval [T n ,T n+ i]. 

The sequence T n is strictly increasing. If it is bounded, then the latter relation above 
yields = 1/2. Hence, T n —> 00 when n — > 00, completing the proof. □ 
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5.2 Stability 



Proof of Theorem 2.2±. Here, we apply [251 Corollary 2.8] to compare the solutions p\ and 
P2 of (|2,4p . Let k be as in (|5.4p . As in the proof of Theorem 12.31 (v)-(i/)~(I) ensure that the 
set of hypotheses [25, (H1*)-(H2*)-(H3*)] are satisfied. Thus, for any t ^ 0, we obtain 

IM*)-P2(i)|| Ll 

^ ||P0,1 — P0,2|| L i 



e /c 'TV(p ,i) + NW N \\ qi \\^ {[0M / e fc ('- r ) / || V div (v x (x) + l x ( Pl )) \\ dx dr 



o 



+ 




'0 

< ||po,i - P0,2|| Ll + 



q[(p) (vi(x) +Xi(pi)) - q' 2 (p) (u 2 (x) + l 2 (p2)) 



qi(p)div (v\{x) +Xi(p 1 )) - q 2 {p) div {u 2 {x) + X 2 (p 2 )) 



dr 



dx dr 



L°°([0,fl]) 

TV (po,i)+tNW N \\ qi \\ Lao{[0:R]) (||Vdiv^|| Ll +C/(||po,i|| L i) 



*lki ~~ ^2 1 1 L o 
+t lk2|| L o 



([0.JJJ) 



([0.JJJ) 



+ P2|| L °°([0,fl]) 
+*IMIl°°([0,.R]) 



+Ji(p 1 )|| Loo 



[ ^(pi(r)) -X 2 (p 2 (r)) 
•/ o 



dr 



div (i/j - z^ 2 +Xi(pi) -X 2 (/9i)) 
- 92|| L -.([0,iJ])ll div (^l + Z l(pl))|| L i + ll®llL»([0,fl3) 

Hence, denoting 
a(t) 
bit) 
7i 

72 



div (X 2 (/9i) -X 2 (p 2 )) 



dr. 



TV(po,i)+*AW7v|| ftll^M) M| V div z^i || L i + C/(||po,i|| L i" 
^ (ll^llL-do^])^* ^) + IMIl~([o,r])) 

h'i - ^||l»([o,hi) (h^IIl- + ^MIIauIIlO) + IKIL^o^lh - 

+ ll^llL»fro^ sup \\ x ^p) - X 2(p)|| L oc 



lift ~ 92 IIl°°([0,«]) 

||div^|| Ll +Cj(|| P0,l|| L i 

+ ll?2||L«([0,fl]) 



div(i/j — 1/2) | L i + sup |divli(p) - divX 2 (p)|| L1 

pGL^R^O.ii]) / 

we get 

||pi - P2(*)|| L i < ||po,i - Po,2|| L i + l\te kt a{t) + tj 2 + b(t) / \\pi - p 2 (r)\\ dr . 

J 

By Gronwall Lemma, we obtain 
\\f>i - P2(*)|| L i < ||po,i -Po,2|| L i +jite kt a(t) +t7 2 

+b{t)e^ du f e~ X b(u)du (\\po,i -Pod Ll +7ire fcr a(r) +r 72 )dr 
J 
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Noting that 



we get 



Jo b ^ du + b(t) f e~ Jo" b (") du dr 



m 

6(0) 



||pi-P2(i)|| L i < (||P0,i-P0, 2 || L1 +7iie H a(i)+i72) ^|ye* fe W 
completing the proof. □ 



5.3 Driving examples 

Proof of Proposition \&.H Let p solve (jl.ip in the sense of Definition 12. 11 Define w(t, x) = 
v{x) + (l{p(t))\ (a?). Then, p e C°([0, T];L 1 (R Ar ; [0, i?])) is also a weak solution to (I2TD . 
Then, set as above = pv(p) and, for any <p € C£°(R x R^jR), we have 



(/? dt<p + q(p) w(t, x) ■ V<p(t, x)) dx dt = . 

Fix positive T and M. For a > define Y Q E C°°(R; R) so that Y a (s) = for s ^ 0, Y a (s) = 1 
for O a and Y^(s) > for all s. Let = Y a (t -a)- Y a (t - T). 

Call K M = B(0,M) \ Set p £ C~([-l, 1]; [0, 1]) such that f R p(r) dr = 1. For <& > 0, 
let p$(x) = ^ p, (\\x\\ and define -0 = * the convolution being over all of R^. 

Note that = on R N \ B(dK M ,0). If a? € B(dK M ,tf), we have 

VV>(aO = / V^?(x - y) dy = - / ^(x - y) h{y) dcr(y) , 
•/-fir M ./afsr M 

where n(x) is the outer normal to Km at any x £ 8Km- Then, we choose </?(£, x) = £(t) VK^)- 
Clearly, </? G C£°(R x R^jR). The definition of weak solution gives 

0=/ / (p(t,x)('(t)ip(x) +q(p(t,x)) w(t,x)((t)Vip(x)) dxdt . 
Jul Jr n ^ ' 

Letting a — > 0, we obtain 

= / (p(0, x) — p(T, x)) ^(x) dx + / / q(p(t,xf) w(t,x)\7ip(x) dx dt 
Jr n ' Jo Jr n 



[ (p(0,x)-p(T,x))^(x) 
Jr^ 

(p(t, x)) w(t, x) • h(y) p${x - y) da(y) dx dt 



dx 

7 / 

'0 JR N JdK M 

(p(0, x) — p(T, x)) ip(x) dx 

f [ q(p(t,y + ^z))w(t,y + ^z)-n(y)p(\\z\\)dzda(y)dt. 
JdK M JB{Q,\) 
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That is to say, setting H#(t) = {(y, z) : w(t, y + tfz)- h(y) < 0} and dQ M = dtt n B(0, M) 
(p(0,x) - p(T,x)) ip(x) dx 

Q {p{t, y + ■&%)) w (t> V + ■&%)■ h(y) p{\\z\\) dz da(y) dt 




>0 JdK M JB(0,1) 
t-T 




'0 JdK M J -8(0,1) 



X H , t) {y,z)<l[p{t,y + tiz)) w{t,y + fiz) ■ h(y) fi(\\z\\) dz da{y) 



dt 



Jan M jb(o,i) 



X H ,Ay,z)q{p(t,y + dz)) w(t,y + &z) ■ h(y)p(\\z\\) dz do{y) dt 



T 



+ 



JdK M \dn M -JB(0,l) 



X H ,Av, z) Q {p(t,y + -dz)) w{t,y + $z) ■ h(y) p(\\z\\) dz da{y) dt 




10 Jan M JB(o,i) 



X H (t) (y^ z ) \W\\-L°oW(t,y + $z) ■ h(y)p(\\z\\) dz dcr(y) dt + m 



where m can be arbitrarily small provided M is sufficiently large, since p G L 1 ([0,T] x 
M. N ; [0, R]). Letting 1? — >• 0, M — >• 00 and using (|3.2|) . by the Dominated Convergence Theorem 
we get 

0</ (p(0,x)-p(T,x))dx = - p(T,x)dx, 

since by hypothesis spt po C ^. Besides, the positivity of p gives us — f cn p(t, x) dx ^ 0. 
Finally, we have J cn p(T, x) = and spt p(T) Cfl. □ 

Proo/ of Lemma HOI Let p G L 1 (R Ar ; [0,i?]). Note that 1 can be rewritten as Z{p) = 



—e p* (Vr/) / y 1 + ||p * V?7|| 2 . By the assumptions on r\ and the properties of the convolution 
product, we have 

-l|Vr?|| Loo , 

||Z(P)|| L 1 < e||p|| L i||Vr?|| L i, 



ll L 

||VX(p)|| Loo < e||p|| L i||V»7|| wl ,oo + 

||VZ(p)|| L1 < e||p|| L i||V»7|| wl ,oo (l + R 2 \\Vr,\\ hl 
Hence (1.1) is satisfied. Let us check now (1.2) 



V 2 r? 



L 1 



L 1 



V 2 X(p) li <e J R 2 ||p|| Ll ||V7 ? || w2 , 1 ||Vr ? || 2 ^ 1 , 1 (l + 4 + 3i? 2 ||V7 ? || 2 ^ M ). 



Let n,r 2 e L l (R N ; [0,R]). We have: 



||X(n)-X(r 2 )|| Loc ^ e^l + i^llVT/II^JIlVT/H^llra-rall^, 

||j(ri)-2:(r 2 )|| L1 < ^(l + ^IIVT/ll^llVT/ll^lln-^llLi, 

||div (X(n) -X(r 2 ))|| L1 < e||r 1 -r 2 || Ll ||Vr ? || wl>1 (l + 8i? 2 ||Vr 7 ||^ M + 3^1^7711^), 
completing the proof. □ 
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Proof of Proposition 13.31 The fact that (v) holds follows through simple computations 
from (fi), (g), (17), (ct) and the assumptions on d£l or a,dgQ. Condition (I) follows from 
Lemma 13.21 

To prove the invariance property (P), we verify that 



y g(x) + S(x) + (Z(p)) (x)J • n(x) > 

for all x E d£l. By (g), (|3.5|) and (|1.2p . it is sufficient to prove that 

(p * Vt])(x) ■ n(x) 



A ^ e sup 



Note that 

(p * S7rj)(x) ■ n(x) 



< sup (p* Vt?)(x) • n(x) < ||p|| LOC ||Vr/|| L i ||n|| LOO < -R||Vr?|| L i 

yi + || (p* Vt7)(s)|| xedn 



sup 

xedn 



completing the proof, by Proposition 13.1 



□ 



Proof of Lemma \8-4[ Note that 1 as X{p) = e J RN p(y) V [r](x — y) g ({y — x) ■ v{x))j dy. 
Let r v be such that B(0,r v ) D sptr?. Fix p € L 1 ^; [0,i?]). Let us check (1.1): 



11^)11 



L°° 



lLi|l r ?llwi>-ll5'llwi.- i 1 + (! + r r?)ll^llwi>-) 



^ e IIpIIli|I 7 /IIwi.iNIwi.«' ( 1 + ( 1 + r »j)ll 1/ llwi.«') > 

ll V ^)Hl 



< e 



| L l||//|| W 2, o||y|| W 2,oo 



1 + 2{r v + l)||i/|| w i,co + (1 + r^lMlwL- + ( 2 + r v) II Vz iw%= 



||VX(p)|| rl < e||p|lLl|kllw2,l|bll W 2-°o 



1 + 2(r„ + l)||i/|| w i,oc + (1 + r^lMl^i.oo + (2 + r„)||Vi/|| w i,o 



Passing to (1.2): 



V 2 X(p) 



L 1 



Ili|I t /IIw3.i||0IIw3.°° 



1 + 3(1 + r n )||i/|| w i,oo + 3(1 + r^fPW^i.o* + 3(2 + r v )\\Vu\\ wl , 



+ (1 + ^) 3 ||l / ||w 1 .°° + 3 ( 2 + r r?) 2 H zy llw 2 .°° + ( 3 + r r) 



V 2 i/ 









Let ri,r 2 £ L 1 (M Ar ; [0, R]). The operator X is linear in p, hence 



|Z(ri)-X(r 2 )| 



< e n — r 2 



Li|i'/iiwi.°°iiy|iwi.°° 



(l + (l+r„)||i/| 



||X(n) -X(r 2 )|| LOO «S e||n-r 2 || L i||?7|| w i,i||5|| w i, t x,(l + (1 + r,,)||i/|| w i 1 ao) 



||div(X(n)-X(r 2 ))| 



^ e Iri — r 2 



Lill'/IIW^illi/IIW 2 ^ 



1 + 2(1 + 

IMI W 1 ' 00 "I" (1 "I" r v) IMIw 1 - 
+(2 + r,)||Vi/|| wl>0 
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showing that (1.3) is satisfied. 

Proof of Provosition[S75\. Fix p G L 1 (M JV ; [0,R]). Let us check (1.1) and (1.2): 



11^)11 



IliII 7 /IIwi.«IIvIIwi.'» i 1 + i 1 + i )\\9\\w^) > 



||^)|| L i < eiipii^ikiiwMii^iiw^Ci + a+^Niw^) 

|VX(p)| 



L°° 



lw 2 > c 



x ( 1 + 2(1 + l)||^|| wl)00 + (1 + ^) 2 ||5llwi.~ + (2 + ^)||Vg|| wl . 

i + R 2 MlvMw^ (i + (i + ^)IMIwi.~) : 



|VZ(p)|| L1 < e||p|| L i||»7|| wa ,i||y|lw».' 



x i + 2(* + l)||0|| wliOO + (l + . 



lw^ + (2 + ^)||V 5 || 



1 + 



iWlMlw^ (l + (l + ^)||5ll wl -) 2 ) 



Passing to (1.2): 



L 1 



IlMI'/IIws. 



wi.oo + 3(1 + £) 2 \\g\\ 2 w i, x + 3(2 + ^)|| V 5 || wl „ 



1 + 3(1 + 

+(1 + ^llsllw^ + 3(2 + ^) 2 ||9|lw2, 0O + (3 + t) 

i + R 2 IkllwM IMlV.« (i + (i + ^)ll5ll w ^) 2 



V 2 5 



w 1 - 



+3eJ2|M| w i,i|Ml w v» (! + (! + 

x ll^llLilhllwa.i ll^llwa.- 

x fl + 2(* + l)|M| w i„» + (1 + ^) 2 ||5llwi,oo + (2 + ^||V 5 || wllC 



i + R 2 hllwM Ibllwi- (! + (! + ^)lbll w ^) : 



lLi|l'/llw3.i 



W"3,oo 



1 + 3(1 + £)\\g\\ w ^ + 3(1 + £) 2 \\g\\ 2 w ^ + 3(2 + £)\\Vg\\ wl „ 



+(1 + ^) 3 ||5llwi- + 3(2 + t) 2 \\g\\ 2 w s,oo + (3 + £) 



V 2 5 



w 1 ^ 



+31?|M| w i, 1 |M| w i, o (l + (l+^)l^llwi.-) 

1 + 2{£ + l)||^|| wl>00 + (1 + ^) 2 ||5llwi,oo + (2 + £)\\Vg\\ w i„ 



i + R 2 IkllwM IMIw^ (i + (i + ^)ll5ll w ^) : 



□ 



Finally, (1.3) is proved exactly as in Lemma 13.21 

To prove the invariance property (P), following the same procedure as in the proof of 

Proposition 13.31 we check that (u(x) + (l(p)) (x^j ■ n(x) ^ for all x G d£l. In fact, if 
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x g <9ft, 



v(x) + (x)\ ■ n(x 



I 



g(x) + A n(x) — e 



\ 



YJq p(v) ~ y) ^{(y - x ) • d ^ 

^1+" V/^p(y)r?(x-y)v?((y-j;) -5r(x))dy 



nix) 



dy 



> \(x) + g(x) • n(x) - s WcpWi^ p(y) V (rj(x - y)) ■ n(x) 

Jn 

-e / P(y) »7(z - y) ((y - x ) • a( x )) [(V ~ x ) n(x) - 5(2;) • n(a;)) dy 

./n 

^ A(x) -eE||^|| LO0 ||Vry|| L i -ei^y/IL^ / n(x - y) \(y - x) Vg(x) n(x)\ dy . 

Jn 

> A(x) -eR(\\ip\\ Loa ||V??|| L1 + \\<p'\\ LOO \\ri\\ L i £ ||V 5 || Loo ) 
completing the proof, by Proposition 13.11 

Proof of Proposition [WTfA Denote I = diamsptr/. Straightforward computations give: 



□ 



II^WIL- < £ 

|VX(p)|| Loo < £ 

II^ILi < £ 

||VZ(p)|| L1 < £ 

V 2 l(p) <: : 



pIIli II V7 ?IIl°° IMIl« i 

V 2 ?7 



^IIl°° + ll Vr ?llL- IML« (IIsIIl* +^11^1^°°)) 



pIIli II v? ?IIli IMIl°° » 



pIIli 



x llvlltoo + l|v??|| L i ||^'|| LOO (Hs-IIl- +^ll v 3ll L -)) 



pIIli l|V77|| W 2,i|M| W 2,oo ( 1 + 3||g|| Loc + ||V5r|| Loo (3^ + 2) + 
giving (1.1) and (1.2). The proof of (1.3) is immediate by the linearity of I in p. 



V 2 , 



□ 



A Appendix: Geometrical Issues Related to Q 

The framework presented in Section [3]can be adapted to various real situations. For instance, 
the "walls" <9ft may stop the visibility of the individuals. Then, from a modeling point of 
view, it can be reasonable to introduce the set 



{y G ft: x + a(y - x) G ft , Vcr G [0, 1]} 



(A.l) 



of the points in ft visible from x. Correspondingly, the nonlocal operator X in (|1.2p can be 
modified intending the convolution as follows: 



(p * n) (x) = I p(y) n(x- y) dy 



(A.2) 



The above relation means that the individual at x evaluates an average of the densities p{y) at 
all values y that are visible from x. With these choices, the validity of condition (I) essentially 
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depends on the geometry of f2. In particular, if is convex, then Q, x = Q and (I) holds by 
Proposition 13.31 

Here we only show how to choose a discomfort so that (P) holds for (|l.ip - ([L2"|) - (|A.2| ), 
The case of the nonlocal operator (|1.3|) is entirely similar. 
To this aim, introduce the set 

H x = {y G n : Vr?(x - y) ■ n{x) ^ 0} . (A.3) 

Proposition A.l. Assume that (CI) and (rj) hold. Let I be as in where the convolution 

is intended as in (|A.2p . // 

5(x)=sR Vr](x-y)dy for x G dQ (A.4) 

then (5(x) + l(p)(x)) ■n(x) ^ for all x G d£l and p G L^M^; [0, R]), so that (P) is satisfied. 

If moreover r/(x) = f](\\x\\) with f) G C^([0, +oo[ ; R + ) and f)' ^ 0, then £l x C H x along dil 
and the discomfort can be defined on all Q by 

5(x) = sR Vt/(x — y) dy 
Jn x 

satisfying to property (P). 



Proof. We have l(p)(x) = —e f n p(y)Vi](x — y)dy j y 1 + ||Vp*7/|| . Hence, to ensure 
that (5(x) +l{p)(x)) ■ n(x) ^ for all x G d£l and p G L 1 ^; [0, R}) we require 



5(x) ■ n(x) ^ e sup / p(y) Vr/(x — y) ■ n{x) dy . 

peL 1 (R JV ;[0,ii])^n :c 

In the latter expression, the supremum is attained for p = R \^ nH , since £l x n H x = {y G 
&x ; Vr/(x — y) ■ n(x) ^ 0}. It is thus sufficient to define 8 so that 

5(x) ■ n(x) ^ e R Vr](x — y) ■ n{x) dy , 

J Q x nH x 

proving that (|A.4p implies (P), by Proposition 13.11 

If r](x) = f](\\x\\), then Vr/(x) = ?/(||x||) -Su. Since fj' ^ 0, 

H x ^ {y(=n x : (x-y)-n(x) ^ 0} 

= jy G O x : y = x + an(x) + uj with a ^ and u G r^x) -1 j 

The latter inclusion holds since Sl x is a convex set contained in Q and containing x. □ 
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The choice (|l,2p - ()A.2|) is appealing from the modeling point of view, but not easily 
tractable from both the numerical and the analytical points of view, without major restrictions 
on the geometry of fi and on rj. Therefore, we consider also the following choice. 

Proposition A. 2. Assume that (ft) and (r\) hold. Let X be as in with the usual 

convolution ()3.3|) . // 



then (5(x) +I{p){x)) ■ n(x) ^ for all x € 89. and p € L 1 ^; [0,R]), so that (P) holds. 
With this choice, a way to define 5 on all Q. could be 



where a is as in (a). 

The proof is entirely similar to that of Proposition IA.1I 
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